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We discuss the supersymmetric formulation of the nonhermitian (3 = 2 random matrix partition 
function with one bosonic flavor. This partition function is regularized by adding one conjugate 
boson and fermion each. A supersymmetric nonlinear cr-model for the resulting Goldstone degrees 
of freedom is obtained using symmetry arguments only. For a Gaussian probability distribution the 
same results are derived using superbosonization and the complex orthogonal polynomial method. 
The symmetry arguments apply to any model with the same symmetries and a mass gap, and 
demonstrate the universality of the nonlinear cr-model. 

o 

^ ' I. INTRODUCTION 

G ' 

, There exists a vast literature [l[ showing that the spectra of many physical systems on the scale of the average level 
spacing (the microscopic scale) are correlated according to universal laws given by random matrix theory. They can be 
classified according to their (anti-)unitary symmetries and invariant bilinear forms 0, E|. The reason for this universal 
behavior is that many physical systems and random matrix models alike can be reduced to field theories with only 
Goldstone degrees of freedom. On general grounds, such a theory is a nonlinear cr-model and is determined uniquely by 
the pattern of symmetry breaking and convergence requirements. One can therefore obtain the microscopic correlation 
functions from symmetry arguments alone without ever referring to a random matrix model. The construction of the 
nonlinear tr-models from symmetries is standard and well known for Hermitian systems [U, HI • 

For nonhermitian problems, which appear in open quantum systems, e.g., in the theory of 5-matrix fluctuations, or 
Euclidean QCD at nonzero chemical potential, the situation has been investigated to a much lesser extent. Here the 
formulation of cr-modcls based on symmetries has been studied in detail only for the case of partition functions given 
\ by products of determinants (i.e., fermionic theories) [Ei fl 0, Hi • A major difference between fermionic theories and 
bosonic or supersymmetric ones is that there are no convergence problems for Grassmann integrals. For bosonic non- 
' hermitian partition functions (with inverse determinants) only very few derivations of a cr-model from the underlying 
\^ , symmetries can be found in the literature; in fact, the only works known to us are (9l. [Io|. In the supersymmetric 
£SJ ' case (i.e., with both fermions and bosons) we are aware of only one model - the generating function of the spectral 
, density of a Hermitian ensemble deformed by an antihermitian ensemble - where a cr-model for the Goldstone degrees 
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of freedom [HI, [13, [3 has been obtained. However, that cr-model was derived by direct calculation, not by symmetry 
arguments. The main objective of the present paper is to show that a cr-model with exclusively Goldstone degrees of 
<^> freedom can be obtained from symmetry arguments alone also in supersymmetric cases. 

In this paper we study the symmetry class whose simplest representative is a model where a complex Hermitian 
Gaussian random matrix ensemble is deformed by a complex antihermitian Gaussian random matrix ensemble. This 
' model was introduced by Fyodorov, Khoruzhenko and Sommers [ll| and for this reason it will be called the FKS 
model. Its introduction was motivated by the study of the distribution of resonance poles for systems with broken time- 



reversal invariance [141 ]. More recently it was used to describe the Hatano-Nelson model [13| and, in its unquenched 
form, QCD in three dimensions at nonzero chemical potential 0, [H[ • 

In this paper we study the FKS model [ll[ for one bosonic flavor, the partition function of which is defined by 

z-^;a) = ( -—- i- , A (i) 



, Det(z + H + A) 

The average is over the Gaussian probability distribution with distribution function 

P(H,A)= C -^ H2 -^ A \ (2) 

with H a Hermitian N x N matrix and A an antihermitian N x N matrix. Note that P(H, A) is invariant by a 
unitary change of basis, H \— > gHg^ 1 , A \— > gAg^ 1 , g e U(iV). As was argued in [l(j, there is a major difference 
between fermionic and bosonic partition functions: in the large- N microscopic limit, the partition function with one 
fermionic flavor does not de pen d on the nonhcrmiticity parameter, whereas the bosonic partition function |T]) does. 
This behavior was found in [1 01 ] using the method of complex orthogonal polynomials for a random matrix model of 
QCD at nonzero chemical potential [1 61 ] - Its explanation was based on the observation that the partition function ([1]) 
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needs to be regularized. In the context of a cr-model formulation, the technical reason is that an inverse determinant 
of a nonhermitian matrix cannot be written as a Gaussian integral in general. The regularized partition function (To| 

Z- 1 (z}\z,z*;a) = ^(Det(z} + H-A)Vet- 1 (* + i 2_ A 2 + ^ + ^)) (3) 

reduces to the partition function for one bosonic flavor for z*j — > z* . It has flavors with opposite charges resulting in a 
ground state which rotates as a function of the nonhermiticity parameter. For one fermionic flavor no regularization 
is necessary and the ground state does not rotate, so that the free energy does not depend on the nonhermiticity 
parameter. The regularization procedure of the inverse determinant is known as Hermitization [l3l ITtI fl8j ■ 

The partition function ([3]) is well understood for zj — > oo in which case it is the two-flavor phase quenched bosonic 
partition function. In that case, because of a complex conjugated singularity, it diverges logarithmically with e 
In cr-model language, the singularity is due to a Goldstone boson with a mass that vanishes as e for e — * P3. Also, 
the partition function ([3]) for e — > acquires a Goldstone fermion of mass z* — zj . We therefore expect the behavior 

Z_ 1 {z}\z,z*;a)~(z*-z})\oge + O(e ). (4) 

In the present paper we will derive this result from the cr-model for the microscopic limit of ([3]). We will also show that 
the C(e°) term agrees with the partition function ([T]) evaluated by the method of complex orthogonal polynomials. 

A major issue with cr-models for nonhermitian random matrix models is the proper choice of integration manifold. 
There are two important developments which have made this choice much less ad hoc. The first of these was the 
introduction of the so-called Ingham-Siegel integral [J] as an alternative of the Hubbard- Stratonovich transforma- 
tion. This work provided a simple explanation of the structure of the integration manifold for inverse determinants 
when convergence arguments are essential. The second development was the introduction of superbosonization which 
extends the Ingham-Sie gel ap proach to include fermions in a unified fashion [20I.I21I]. Earlier versions of superbosoniza- 
tion appeared in [22|, HH, H3 HE HH , but the method was put on a mathematically rigorous footin g on ly in [20I l2l| . 
Although it is straightforward to obtain the correct a- model using the superbosonization formula of [2 (J . |21[ , it can be 
a technically challenging task to evaluate the integrals for more than a few degrees of freedom. The superbosonization 
method was applied to nonhermitian chiral random matrix ensembles in [27l |. However, the construction of a nonlinear 
cr-model containing exclusively Goldstone degrees of freedom was not given in that paper. 

In this paper we will derive universal results for the symmetry class of the FKS model relying on symmetry 
arguments only. This makes it manifest that our results apply to all models in the same symmetry class and a mass 
gap. Analytical results for finite N cannot be obtained from general arguments and require a detailed calculation. 
Such results will be derived for the FKS model with Gaussian probability distributions using two independent methods, 
the superbosonization method and the complex orthogonal polynomial method. Each method has its own merits and 
both deserve to be discussed. In particular, relations between partition functions with different degrees of freedom 
appear naturally in the complex orthogonal polynomial method. Using superbosonization the universal cr-model can 
be recovered from the finite- N results by taking the microscopic limit and eliminating the massive modes. 

Finally, let us mention that the quenched spectral density of the FKS model has also been derived @ by means of 
the replica limit of the Toda lattice equations [28] . It was shown that the quenched spectral density is the product of 
a fermionic and a bosonic partition function. To derive the spectral density for one fermionic flavor using the Toda 
lattice hierarchy, one needs precisely the partition function ((3]) [IH . 

In this paper we will first derive the cr-model for the microscopic limit of (J3]) using symmetry arguments only 
(Section II). Results for finite N will be derived using the superbosonization formula (Section III) and the complex 
complex orthogonal polynomial method (Section V). In Section IV we will recover the universal cr-model from the 
finite- V results for the FKS model. Concluding remarks are made in Section VI. 



II. SYMMETRIES AND SUPERSYMMETRIC u-MODEL 



In this section we will derive the universal supersymmctric cr-model for the symmetry class of which the partition 
function ([3J is the simplest representative. The derivation is based on the symmetries of ([3]) only, and is valid for any 
other model in the same symmetry class. We will first consider the bosonic sector with one bosonic flavor 4> + and one 
conjugate bosonic flavor 0_ . An earlier study of this sector was made in Q. 
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A. The Phase Quenched Bosonic Partition Function 



The regularized phase quenched bosonic partition function can be written as 

Z pq _ bos (z, z*; a) = ^Det" 1 ( z « + % _ A * + l + A )) ■ (5) 

We will evaluate this partition function in the microscopic limit, keeping TVImz and Na 2 fixed as N — > oo. To start 
the argument, we cast the inverse determinant in the form of a Gaussian integral: 



2V 



k=l 



. hos (z, z*; a) = I J] dftdtlW-WJ (exp i ( <t>* + <t>*- ) I z , + % A * + \ + A ) ( ) ) • (6) 



After averaging over H and A the partition function can be expressed as an integral over the positive Hermitian 2x2 
matrix Q of U(./V)-invariant variables 



JV 



k=l 



(7) 



If we were dealing with fcrmions and compact symmetries, we could now consider a form of 'maximum flavor symmetry' 
using the theoretical arguments of Peskin 0, Q ■ A nonzero expectation value of such a form would signal spontaneous 
symmetry breaking in the thermodynamic limit. Our analysis must be somewhat different, however, as we are facing 
the case of noncompact bosons. To express the microscopic limit of the partition function (J5J) in terms of the Goldstone 
degrees of freedom residing in Q, we will require that the resulting integration measure and Lagrangian have the same 
transformation properties as the corresponding objects of the original partition function. 



The complex group GL(2) acts on the matrix Q and the vector variables (j) = ((/)+ 4>-) an d 4>* — \ J a s 

r^94>*, <^0<? + , Q^gQg f , g e gl(2) . (8) 

Under such transformations the integration measure Jlfcli d<jfid<j)^d<fi+d(j)*]' gets multiplied by the Jacobi determinant 
|Det g\ 2N ■ Thus, arguing by symmetry and equating the transformation behaviors, the corresponding measure in the 
©/-variables is inferred to be Det N (Q) dQ where dQ denotes a GL(2)-invariant measure for Q. Note also that the 
transformation law Q t— > gQg^ preserves the properties of Hermiticity and positivity of the matrix Q. 

For z = z*, a — , and e^Owe see that the partition function §E§ is invariant under the subgroup G C GL(2) of 
flavor transformations Q i— > TQT^ which preserve the Hermitian quadratic form 

<t>+<t>- + 010+ = TrQcri , iti = (J J) ' ( 9 ) 

of the bosonic degrees of freedom <p . Equivalently, the matrices T £ G are subject to the condition T^o\T = o\ . Since 
the Hermitian quadratic form determined by a± is of signature (1, 1) — to see this, one makes a unitary conjugation 
transforming o\ into 0-3 = diag(l, —I) — the symmetry group G of our problem is identified as G = U(l, I). 

Next we search for the manifold of Goldstone degrees of freedom (or the target space of the nonlinear er-modcl) 
inside the space of matrices Q. For that purpose, we make a temporary change of variables from Q to X := \Qo\ . The 
new matrices X satisfy X = —u\X^a\ and thus lie in the Lie algebra Lie(G) of G = U(l, 1). Now, using Q 1— > TQT^ 
and the relation T^a\ = (TiT^ 1 we see that G acts on X £ Lie(G) by the adjoint representation X 1— > TXT^ 1 . 

The space of 2 x 2 matrices Q subject to the conditions Q = > is a cone of real dimension four. Writing 

t:t~) =Q] ' < io) 

this cone is given by the inequalities Q++ > 0, Q > 0, and Q |_ = \Q^ | 2 < Q++Q . Thus our new 

matrices X = iQ<Ji do not occupy the entire Lie algebra of G but lie in the cone C + C Lie(G) which arises as the 
corresponding image by the map Q ^ iQui . The cone C + may be viewed as the 'space of states' of our problem. 

Let us look at the positive cone C + in a little bit of detail. First, notice that Lie(G) is generated as a Lie algebra 
over the real numbers by the four generators il, iui , a 2 ■, an d 03 . Let K C G be the maximal compact subgroup which 
is generated by the first two, il and ici . Thus K is the group U(l) x U(l) of elements e^ a+/3cri ' with a,/3 £ [0,27r]. 
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Now consider t + := C + n Lie(if), the intersection of the cone G + with the Lie algebra of K . One may ask whether 
the elements £ £ Lie(G) can be conjugated into ? + by the adjoint action £ i— > T£T _1 of G. The answer is that this 
is not possible in general, because Lie(G) is the Lie algebra of a noncompact group. Nevertheless, the cone G + does 
have the special property that each of its elements is conjugate to some A £ t + by the adjoint action of G. (This 
follows from basic principles of linear algebra and Lie theory and can be easily verified by direct calculation for our 
simple case of 2 x 2 matrices.) Thus each £ £ C + can be presented in the 'diagonalized' form 

£ = TXT- 1 , T £ G , A G 6+ . (11) 

To introduce the proper mathematical language, we say that each orbit of the adjoint G-action on G + hits the slice 
6+ C Lie(K) at least once (actually, exactly once). It follows from the diagonalization ifTTj) and the abelian nature 
of K = U(l) x U(l) that C + has the structure of a direct product (G/K) x t + . (More generally, in the case of a 
nonabelian group K, the cone C + is an associated bundle G Xk £ + -) Moreover, since K C G is a maximal compact 
subgroup, the quotient G/K is a symmetric space of noncompact type. In the case at hand we get the identification 

G/A" = U(l,l)/U(l)xU(l) = H 2 (12) 

with a two-dimensional hyperboloid H 2 . To summarize the present discussion: the space of states of our problem, the 
positive cone G + , has a decomposition (mathematically speaking, a 'fibration') G + = (G/K) x t + by adjoint G-orbits 
all of which are isomorphic to the same noncompact symmetric space G/K . 

Now, by some dynamical principle beyond the reach of symmetry arguments, the system selects one of the G-orbits 
of the fibration C + = (G/K) x t + for its Goldstone manifold or vacuum orbit. This G-orbit will in general be specified 
by an clement A £ t + (g> C of the complcxification of t + . In the case under consideration we have 

A = iA l + iAio-i , (13) 

where Ao and Ai would have to be real numbers (with |A | < Ai) in order for A to be in t + , but in view of the principle 
of steepest descent (or deformation of the integration contour into the complex plane) we should be prepared for Ao 
and/or Ai to deviate from the real axis. For example, in the case of the Gaussian ensemble ([2]) one finds that the 
Q-integral for z = , a = , and e-*0 has a saddle point (a maximum of the integrand) at Ao = and Ai = 1 , with 
Ao becoming imaginary as Re z moves away from zero. In general, Ao and Ai take some other values. While these 
may be hard to compute, we will see that the universal results emerging in the microscopic limit do not depend on 
them. 

Returning to our original notation, we have identified a Goldstone or saddle-point manifold of matrices Q: 

Q = TQ T^ , Q = X 1 l + X a 1 , T £ G . (14) 

From the discussion above, we know that this G-orbit Q = TQqT^ is always isomorphic to the quotient G/K of the 
noncompact group G by a maximal compact subgroup K C G. 

There exist very many ways of parameterizing the G-orbit Q = TQqT^ . One possible choice is by a diffcomorphism 
H 2 ~ R 2 , exponentiating the real plane R 2 spanned by the generators <72 and 173 as follows: 

T = e ua 3 /2 e s<T 2 /2 ^ q = T Q oT t = XiC ua 3 /2 c sa 2c ua 3 /2 + ^ > ^ = /0 -l\ ^ g g ^ ^ 

Let us briefly pause to mention the following heuristic confirming the present scenario. Suppose we were to go 
beyond the microscopic limit and construct a nonlinear a-model of spatially fluctuating Goldstone modes with target 
space G/K. To give a sensible definition of the functional integral of such a field theory, we need the target space to 
be Ricmannian. Now, for the case of a semisimple noncompact Lie group G it is a fact of differential geometry that 
there is only one way to get a Ricmannian manifold with a G-invariant geometry: divide G by a maximal compact 
subgroup K. In contrast, the situation for fermions with compact symmetries is very different. There, the fibration of 
the state space by orbits of the symmetry group typically contains orbits of several types, corresponding to a variety 
of nonisomorphic compact Ricmannian symmetric spaces. In that situation, unlike what we are facing here, one has 
to appeal to a postulate of 'maximum flavor symmetry' Q to select the proper type of vacuum orbit. 

We are now getting ready to switch on the perturbations a , e , and z ^ z* breaking G-symmetry. Using the 
transformation law Q 1— ► TQ T^ for the Goldstone degrees of freedom, the partition function ([6]) in the presence of 
the symmetry-breaking terms remains unchanged if we simultaneously transform 



C^T t_1 CT _1 , A^T^^AT- 1 , 



(16) 
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with 



Of course the low-energy limit of the partition function © must inherit the invariance under the transformation fT| 

Here, to proceed, we make the assumption that the low-energy measure which is induced on the G-orbit G/K of 
the global mode (or zero mode) converges to the G-invariant measure, dfi(Q), when the regularization parameter e is 
taken to zero. In the case of a compact symmetry group G this assumption always holds true. However, in the present 
case of a noncompact symmetry G (more precisely: a 'nonamenable' symmetry G, see p9|). the G-invariant measure 
is intrinsically unstable with respect to interactions of the Goldstone modes. This circumstance causes a breakdown 
[29l ] of the standard scenario of spontaneous symmetry breaking, zero mode approximation, and universality. 

Nevertheless, in the microscopic limit, i.e., for weakly interacting Goldstone modes in a small enough volume, the 
said assumption does hold true, and the integration measure on G/K in the limit of a = , hnz = , and e — > , 
is the G-invariant measure dfi(Q). The G-invariance of the low-energy partition function then forces the low-energy 
Lagrangian to be G-invariant as well. For the mass term there exists only a single invariant to lowest order in £ : 

Tr(CQ) = AiTrtCTTt) . (18) 

After averaging, there are no terms linear in a . To order 0{a 2 ) there are two possible invariants: 

Tr(AQAQ) and Tr(AQ) Tr(AQ) . (19) 

While these invariants are independent in general, it so happens in the present case of a single flavor that they are 
accidentally the same. To verify this, one may exploit the paramctrization (fT5j) and the relation e u<T3 a 2 = <J 2 c~ ll<T3 to 
find the expressions 

Tr(o- 2 <3) = 2Aisinhs , Tr(o- 2 Q o 2 Q) = 2A 2 (1 + 2 sinh 2 s) - 2A 2 , , (20) 

which show that Tr(fi2Q 02 Q) — Tr(<72Q) Tr(c2<3) ls a constant independent of u and s. 

Thus we need only include the first invariant in the expression for the partition function. Terms of higher order in 
£ and a 2 do not contribute in the microscopic limit and will not be considered here. We also see that the unknown 
parameter Ao just adds to the low-energy Lagrangian an inessential constant, which will not be considered any further 
here (i.e., we set Ao = 0). The remaining unknown Ai is determined by the eigenvalue density of the system, and we may 
take its value to be Ai = 1 by an appropriate choice of units. Note also that Dct N (Q) = Det N (TT^) = IDet T\ 2N = 1. 
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We thus find that the microscopic limit of the phase quenched bosonic partition function is given by 

Z pq - hos (z,z*;a) = [ MQ) j™CQ-h»*' ^«?<^) , (21) 



which is an integral over the coset space G/K of matrices Q = TT^ with G-invariant measure dn(Q). 

Now the two-hyperboloid G/K = U(l, 1)/U(1) x U(l) is the simplest member of a certain family - the Hermitian 
symmetric spaces - with many wonderful properties. In particular, Hermitian symmetric spaces arc Kahlcr manifolds 
and come with a G-invariant, closed and non-degenerate two-form, uj (the Kahler form). In the case at hand, 

ui = — iTr((TiT _1 dT A T _1 dT) , (22) 

which is clearly invariant under left translations T 1— > gT corresponding to the G-action Q 1— > gQ , and also pushes 
down to a well-defined form on the quotient G/K. Using d 2 = and rf(T _1 ) = —T~ 1 (dT)T~ l the expression for ui 
simplifies to 

< J j = idTt(a 1 T- 1 dT). (23) 

We will shortly use this formula to compute the expression of our G-invariant measure dfi(Q) in suitable coordinates. 

To calculate the integral (j2Tj) we use the paramctrization ([15]), and we note that Tr((Q) = ieTrQ + Im(z)Tr a 2 Q- 
From (|20|) we already have the expressions for Tra 2 Q and Tr(a 2 Q) 2 , and for the remaining term in the exponent we 
find eTrQ = 2e cosh u cosh s . To express the measure dfi(Q) of integration we insert the parametrization (|15p for T 
into (f2"3"|) to obtain 

cj = idTr{a 1 T~ 1 dT) = \ dTr(ai e~ s<T2 <7 3 ) A du = rf(sinhs) A du . (24) 



6 



From this result we can say immediately how the measure d/i(Q) looks in the present coordinates. Indeed, since the 
form uj is G-invariant, so is the integration measure d(smhs)du corresponding to it. Because the measure dfi(Q) is 
determined uniquely (up to multiplication by a constant) by G-invariance, we conclude that dfi(Q) oc rf(sinhs) du. 
Assembling terms, the phase quenched partition function (|2ip becomes 

Z pq _ bos (z,z*;a) = lim / e 2UVIm2 sinh s -Wa 2 (i+ 2s inh 2 s ) ( f e -2W £ co S h« cosh^A d ( sinhs )_ (25) 
£ ^°+ Jr \Jr J 

The inner integral over u diverges as | loge| for e — > 0. The outer integral over s is then a Gaussian integral in sinhs 
which is easily done by completing the square. Thus our final result [9j for the partition function is 

^ pq -bo S ^^;a) = |log e |y^c-^ 2 -f Ira2 ^) ( e -0). (26) 



B. The Partition Function for one Boson 



In this subsection we analyze the partition function To that end, we express the determinant in the numerator 
of ([3]) as a Gaussian integral over a U(iV) fundamental vector ip of Grassmann variables ip k . We then combine ip with 
the boson flavors <p± to form a supervector (f> = (0+ c/>_ ip) with adjoint 

**=U*j. (27) 

The low-energy effective degrees of freedom will emerge from a supermatrix Q of U(iV)-invariants, 

Q = $*$ = ^ + 0* 0_ ^ (fa/> = J2 k etc.). (28) 

Note that the boson-boson block of Q is Hermitian and positive as before. The matrix entry Qff = ^ip of the 
fcrmion-fermion block will acquire a nonzero vacuum expectation value and is treated hence as a complex number. 

Guided by the symmetries of the microscopic theory, we are now going to identify the low-energy degrees of freedom 
and the structure of the low-energy Lagrangian. On general field-theoretic grounds, we expect the low-energy theory 
to be a nonlinear cr-model of interacting Goldstonc modes where the target manifold is a symmetric space. 

To see why the target space has to be symmetric — we briefly recall the argument here — one may invoke Fricdan's 
work [30] on the renormalization of nonlinear models and cr-models, which shows that the quantum loop corrections 
to the target space metric are given by contractions of the Riemann curvature tensor; the one-loop correction, in 
particular, is given by the Ricci curvature. (These results, while derived in the classical setting, remain valid in the 
supersymmetric context.) Therefore, in a low-energy fixed point theory the Ricci curvature of the target space must 
be proportional to the metric tensor. It follows that the curvature has to have the property of being covariantly 
constant which, in turn, is the condition for a Ricmannian manifold to be a symmetric space. This result, which is 
fundamental for the renormalization theory of nonlinear models and cr-models, will presently be used. 

Under the most general linear transformation of the supervector 

S^Gta)- 1 , (29) 
the (a priori) superintegration form = JllLi D<& k D<& k * transforms as 

D$D$* 1 ► D<S>D<5>*SDct N (g L ) SDet"^^) . (30) 

To match this transformation behavior, the space of composite variables Q has to be equipped with the Bcrczin 
measure (or superintegration form) DQ SDet Ar (Q) where DQ by definition is invariant under Q 1— > g^Q (gn) 1 ■ 
The Hermitian quadratic form @ is replaced by the boson-fermion mixed form 




^ = STrQEi , Ei= 100 . (31) 
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The symmetry group of this extended Hermitian form is the pseudo-unitary Lie supergroup G = U(l, 1|1) (a close 
variant U(l, 1|2) of which was discussed in detail in [3ll]). 

Given the symmetry group G = U(l,l|l), we now ask again about the fibration of the space of states by G-orbits. 
For that, we temporarily switch from the supermatrices Q to the related supermatrices X = QSi , on which the 
symmetry group G acts by conjugation: 

QSi i-> T(QEi)T~ 1 . (32) 

We know from Section III AI that by this action the number part of every matrix X can be brought to diagonal form. 

Consider first the generic case of 3 x 3 supermatrices X = QSi with three eigenvalues that all differ from one another. 
The orbit of the G-action on such a matrix is a flag supermanifold U(l, 1|1)/U(1) x U(l) x U(l). Such a space is not 
symmetric (indeed, the Ricmannian curvature is not covariantly constant but varies) and by the renormalizability 
criterion reviewed above, it can be ruled out as a candidate for the Goldstone manifold of vacuum states. 

There exists, however, the possibility for another type of G-orbit, which is realized when the fermion-fermion part 
of X becomes degenerate with an eigenvalue of the boson-boson part. Supermatrices X on such orbits are of the form 



X = T A x A T~\ T e G and |A | < Ai . (33) 




Thus the degeneration occurs in one of two different ways: the boson-boson part of X has eigenvalues Ao ± Ai and the 
vacuum expectation value of Xff = ipi/j may hit cither one of these. In both cases our generic G-orbit degenerates to 

G/K = U(l,l|l)/U(l)xU(l|l), (34) 

where K C G is defined for (V»V) = A + Ai by the equation fcSi&i -1 = Si , and for ('(pip) = A — Ai by 

/0 1 0\ 

fcSifc" 1 = Si , Si = 1 . (35) 

V° -v 

The quotient G/K is a symmetric superspace, and thus satisfies the renormalizability criterion, in both cases. There 
exists no dynamical or other reason (not in the microscopic limit anyway) why one of the two G-orbits corresponding 
to the two vacuum expectation values ('tpip) = Xo ± Ai should be preferred over the other. We also note that these two 
G-orbits are disjoint. The low-energy theory is therefore expected to be a nonlinear cr-modcl with a two-component 
target space, i.e., with one connected component for each of the two vevs. 

Without loss, we now simplify the discussion by setting A = and Ai = 1 as before. The low-energy degrees of 
freedom are then represented by two supermatrices Q and Q', 

Q = TS 1 T- 1 Si, Q' = TS' 1 T~ 1 S 1 , TeU(l,l|l). (36) 

Both Q and Q' run through a symmetric superspace G/K = U(l, 1|1)/U(1) x U(l|l), which has the property of 
being Hermitian. This fact will be of great help in expressing the G-invariant integration measures DQ and DQ' in 
coordinates. Here we just note that SDet N (Q) = SDct Ar (TS 1 T" 1 S 1 ) = 1 and 

SDet w (Q') = SDet^TEiT^Ei) = {-l) N . (37) 

To write the result for the partition function in a concise manner, we introduce a superscript a = ±1 and let 
Q a = Q for a = +1 and Q a = Q' for a = — 1, and we denote the G-invariant superintegration form by Dfi(Q a ). The 
invariance arguments of the previous section still apply. In the microscopic limit we thus find 

Z_ 1 {z}\z,z*-a)= ° N I D^Q a ) C iNsrn ^-^ Na2sTr ^^ s \ (38) 

a=±l 

where the mass matrix £ and the extended Pauli matrix S2 are now given by 

(39) 
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To summarize, we have expressed the microscopic limit of Z_x(Zf\z, z*; a) as an integral over Goldstone degrees of 
freedom only. (Of course, the overall normalization factor cannot be fixed by the arguments in this section.) In the 
next section we will rederive the result (|38|) using superbosonization. 

The calculation of the integral (j3"8")l requires an explicit parametrization of the Goldstone degrees of freedom. For 
this purpose we choose a fcrmion-boson factorizcd generalization of the parametrization (1151) : 



T = T f T b , I>=(J°), V=(j 3 a 1 Y T b =(™i)> W = eWV-/ 2 . (40) 



To motivate this choice, let us observe that 




E 3 := -iSiSa = 0-1 (41) 



commutes with T f . Also, let Ef := Si for a = +1 and := Ei for a = -1, so that Q a Y<i = TEfT" 1 in both 
cases. Note (Ef ) 2 = 1. Looking at Q tT E 2 = iT/T 6 EfT 6 _1 T^ 1 E 3 we see that the expression for STr(<2 CT E 2 ) 2 remains 
unchanged from ([20]) but for the addition of a trivial constant: 

STr(Q' T E 2 ) 2 = 3 + 4 sinh 2 s . (42) 

The mass term has the expression 

STr^Q^ = 2ie cosh s cosh u — azj — i(z — z*) sinhs + a/3((z* — Zj)(isinhs — a) + iee" coshs) . (43) 

To express D^Q' 7 ) in coordinates we recall a few facts from (super-)geometry. If x 1 , . . . , x p and . . . , £ q form a 
system of commuting and anti-commuting local coordinates for a Ricmannian supermanifold M with metric tensor 

g = Aij dx l dx } + B ik dx l d£ k + Qj d^ l dx° + D ki d£ k d£ l , (44) 

where Aij = Aji and Dki = —Dik are even and = Cki are odd functions of the £'s, then the associated superintc- 
gration form is 



n-«*-*f4r-4 f .SD* , ' , (32) 



(45) 



Now our supcrmanifolds M = U(l, 1|1)/U(1) x U(l|l) carry the G-invariant metric 

g = 2STr(T- 1 dT) 2 , (46) 

where the subscript p signifies projection on the tangent space of M at the origin; i.e., A v is the component of A that 
anti-commutes with EJ : 

A p := | {A - E^EJ) . (47) 

Using (|46| in conjunction with (j44l |45|) we can express fl = Dji(Q a ) in any coordinate system of our choice. 

However, in the present situation there exists a better way of doing this calculation. Being a Hermitian symmetric 
superspace, M = M a (for a = ±1) comes with a tensor field J called a complex structure. This means that J 2 = — 1 
and J is an isometry of the metric, i.e., g{Ju,Jv) = g(u,v) for any two tangent vector fields u, v. In our case the 
induced action of J on the one-form (T~ 1 dT) p is given by 

(T- 1 dT) p i ► iE^(T- 1 fiT)p = -i(T- 1 dT) p 'EZ . (48) 

Given the (Kahlcr) metric g and the complex structure J one defines a two-form uj (the Kahlcr form) by the equation 
ui(u, v) = g(u, Jv). In the present case we find 



UJ = UJ 



-iSTr(E^(T- 1 dT) p A (T~ l dT) v ) = idSTT^T^dT) . (49) 



For our purposes, the form a; is a useful object to introduce because uj is easier to express than the metric g and yet 
carries enough information to construct the Berezin measure fl = Dfj,(Q cr ). In fact, if oj is expressed in coordinates as 

w= \ (Aij dx l A dx J + B tk dx l A d£, k + C l3 d£ l A dx j + D kl d£ k A d^j , (50) 
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where Ay = — Aji , Dki — Dik , and = —Cm (due to skewness of the wedge product), then we have 

SDet(^)=SDet(^|) (51) 

as a consequence of the properties of the complex structure J relating the metric g with the Kahler form to. 

We are now in a position to compute D^{Q a ) with ease. Inserting the parametrization (|40|) into (j49|) we obtain 

uj a = i d STr ( 1 dT b ) + idSTr^E^T^T^cff/) . (52) 

The term idSTr(Y,%T~ 1 dT b ) = idTr(a 1 W~ 1 dW) = d(sinhs) A du was already computed in ([2"3]). The new term is 

idSTr^EJT^ 1 ^" 1 ^/) = d(sinh(.s) ad(3 + ia(3da) = (icr + sinhs) da A d/3 + a d(sinhs) A d/3 . (53) 

The last summand makes no contribution to the superdeterminant of the metric tensor (since a 2 = 0) and therefore 
can be dropped for the purpose of constructing Dii(Q cr ). The term proportional to da A d/3 = d(3 A da contributes 
the reciprocal of the analytic square root of — (icr + sinhs) 2 . Thus from Eqs. (|5"Tj) and (|4"B")l we have 

DMQ <T )= 5 SinhS) f dadp. (54) 
1—icr sinh s 

Please be advised that the symbol da in this expression means the derivative da = d/ da , although its meaning in 
the previous equation was that of a differential. Thus we are using the same symbol da for two very different objects. 

III. MICROSCOPIC LIMIT OF THE PARTITION FUNCTION 

In this section we evaluate the partition function (|38p in two different limits. First, we extract the contribution 
that diverges as loge for z*, and second, we compute the result for zj = z* which is regular for e — > . 

A. Contribution of order log e 

Using (|4"3"|) for the mass term and the expression (|4"2"|) for STr(QE 2 ) 2 j the partition function (|3"5)) becomes 

Z^(z}\z,z*;a) ~ r N e~i Na2 V a N e~ m ^f [ rf ( sinh5 )^ dad/3 
1 d J 1 — icr sinh s 

a— ±1 

^ g-2We cosh s cosh u+A^s — 2*) sinh s — 2N a 2 sinh 2 s+iNa/3((z* — Zf) (i sinh s — cr)+iee u cosh s) (55) 

Here we display only the factors of alternating phase; the full overall normalization factor will be inserted below. 
The integral over the Grassmann variables a, (3 yields a factor 

Ne e u cosh s + Nia(z* — Zj)(l — \a sinh s) . (56) 

Since the rest of the integrand is even in u we may replace e" in this expression by cosh u . The w-intcgral with the 
resulting term Ne cosh u cosh s is finite in the limit e — > (see next subsection) . Therefore we may drop this term 
here, as we are after the singular contribution cx loge . The u-integral over the remaining term has the asymptotics 

J due- 2Necoshscoshu = 2 | loge | +O(e ) . (57) 

Doing finally the s-integral by completing the square we obtain the leading term 

Z--i{z* f \z, z*- a) = c N \\oge\ a~V - z* f ) e"^ ta 2 (*/a)-f JVa 2 S m(Nz* f + Ntt/2) + O(e ) . (58) 

The normalization constant cjv is found by keeping track of all constants in the calculations above. Using the formula 
volU(iV)/volU(-/V — 1) = (2tt) n /(N — 1)! and Stirling's approximation for the factorial, we find 

c N = e N ^(N/TT) 1 / 2 . (59) 

This result (f5"5)) will be verified by taking the microscopic limit of the exact finite- N results that will derived in Section 
IIVI by means of supcrbosonization and in Section [VTl using the method of complex orthogonal polynomials. 
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B. Contribution of order e° 

For = z* the logarithmic singularity | loge| vanishes. Indeed, doing the Grassmann integrals over a, (3 and the 
it-integral and then sending e->0wc obtain a finite limit 

lim eNcoshs [ du e u - 2eNcoshscoahu = 1 . (60) 



With the substitution q = sinh s the expression (|55p for the partition function now becomes 

Z- 1 (z*\z,z*;a)~e- S * Na2 V e -*W-*/*) f ^(-O <r-a*«V (61) 

J=±i J 1 ~ laq 



Introducing an auxiliary integration by (1 — \<rq) 1 = J °° c *^ 1<yq ^dt we can do the Gaussian integral over q by 
completing the square. The result of this step is immediately expressed in terms of the complementary error function: 

J 1-iaq V aV2N J 

which is defined by crfc(x) = 1 — erf(x) = (2/y^Tr) e~* dt . Our final result reads 

Z_ 1 (z*|z,z* ;a )=e^( 1 + Q2 >2- 3 / 2 ]T c- iaN ^^CTk(aV2N+^^-] . (63) 

a— ±1 ^ d'S/^N J 

This expression agrees with the result obtained from the microscopic limit of the finite- N results that will be derived 
by means of supcrbosonization in the next section, and by using the Cauchy transform of orthogonal polynomials in 
Section ED 

IV. SUPERBOSONIZATION 

As discussed in the introduction, symmetry arguments alone are not sufficient to derive results for finite N. This 
requires an exact evaluation of the partition function which, in this section, is achieved by the method of supcr- 
bosonization. In the next section we will obtain exact finite- iV results by means of complex ortho gon al polynomials. 

The superbosonization method was introduced to address problems with nongaussian disorder [20l . l2ll , I22I I23I l24l . 
[25l . [26| . Its main idea is to reduce an integral with symmetries to a lower-dimensional integral. To give a simple 
example illustrating this general idea, consider a function / of complex variables z 1 , . . . ,z N . If / depends only on 
x = y^fr— j |z fc | 2 then the integral of / over a U(A^)-invariant domain in C N can be reduced to an integral over just x. 
Similarly, the Grassmann integral of a function /Q^ fe=1 ip k '4> k ) of anti-commuting variables if) 1 * and is known [32j 
to be expressible as an integral of f(y) over y £ U(l). 

Based on results from invariant theory, superbosonization extends this reduction idea to the general case of invariant 
functions of supervectors. In the bosonic sector, the method is equivalent to the one introduced in However, in Q 
the fermionic degrees of freedom were bosonized in the usual way by means of a Hubbard-Stratonovich transformation 
(we will refer to this procedure as the hybrid method), whereas in the superbosonization approach the fermionic and 
bosonic variables are treated on equal footing. From our perspective, a major advantage of the superbosonization 
method is that the integration measure is given by a general formula which can be easily applied to a specific case 
such as the bosonic partition function considered in this paper. In order to execute the integrals, it is essential that 
the paramctrization be chosen judiciously. 

The present partition function has also been worked out in a straightforward way using the hybrid method of [3]. 
That calculation is not more complicated than the supcrbosonization method, but since it does not provide us with 
any additional insights we will not discuss the hybrid method any further. 

The starting point for superbosonization of the regularized FKS partition function ([3]) is the representation of the 
inverse determinant as 

Det -i ( ie ^-ff + ^^j = f ^^'g-^^'^+^^J+i^^.+fly+A^^+i^f/fc+ffH-^O^ (64) 

I Z ~\- H y4. 16 J J 

and the fermion determinant as 

Det(z; + H - A) = J dipdiP e -rt k (**,SM+H kl -A kl W\ ( 65 ) 
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Taking the average over the Gaussian distribution ([2]) of H and A we obtain 

AH kl (r + k <t> l -+<t>*-<i>\-^ k i> l )\ _ e -(l/2iV)STrQ£iQ£i ^ /gg) 



iA, E K0+V_-0*V+ + <AV)\ _ „-(a 2 /2Af)STrQS 2 QS 2 



A 



where Q is the supermatrix Q28[) of U(-/V)-invariant bilincars, and the matrices Si and £2 were defined in (|3 1 [) and 
(|39|) . The quadratic terms in the exponents of ((64)) and (|65|) combine to make up the mass term: 



e -^+ <t>+-e<t>- <t>-+iz<t>+ <t>-+^'<t>- <t>%-iz}^^" _ giSTrCQ _ (gg) 

The method of superbosonization allows us now to introduce the matrix elements of Q directly as the new variables 
of integration. Using a formula proved in [2(| the partition function (after rescaling Q — > NQ) reduces to 

Z-Az*Az,z*-,a) = (2m)- N ™}^\, / DQ SDet iV '(NQ) e iJvslY ^" ^ ^0x^-^/2) STr Q E 2Q s 2 _ (6Q) 

vol U(iV — 1) J 

By the superbosonization step of passing to (|69[) . the precise meaning of Q has been transformed: Q is now the 



supermatrix 



q = (2 bb 2 bf ) 



x++ 


X+- 




X-+ 




£- 


11+ 


V- 





(70) 



where the boson-boson block Q^ = X is a positive Hermitian matrix, Qjj = y £ U(l) is a unitary number, and the 
components of Q/6 = = 77-) and Qbf = ^ = (j^^J are Grassmann variables. The Berezin measure is (20I [2~lj 

DQ = (4n 2 i)- 1 d A Xdy d£ + dq + d£,-dr)- SDet _1 (Q) , d A X cx dX ++ dX— dX+_ dX_ + . (71) 

Following the conventions of [2(| we normalize the flat measure d 4 X so that lim^o 6 2 Jx>o Tl ^ x ~ x ^ d 4 X = 1. 

Note that DQ is scale-invariant; it is also invariant under the transformation QEi 1— > QE'i for E^ given in (135p . 

The expression (|69[) for Z_i is suitable for saddle-point analysis in the limit N — > 00. Since we are considering the 
microscopic limit where and Na 2 are held fixed as N goes to infinity, the symmetry-breaking terms are subleading 
in l/N and can be temporarily neglected for the purpose of finding the saddle-point manifold. If we set £ = and 
a 2 = , the logarithm of the integrand becomes 

JVSTr bgQ - ySTr QEiQEi , (72) 

variation of which gives the saddle-point equation Q~ x — EiQEi . The solutions of this equation form the two disjoint 
supermanifolds Q = TEiT^Ei and Q = TE^T^Ex which were described in (|36[) . Note that the signs of the solution 
in the boson-boson sector are fixed by the condition Qtb > . 

Next we perform the integration over the massive modes in the large- N limit. To handle both saddle-point manifolds 
at once, we recall our notation E^ = £1 for a = +1 and EJ = E^ for a = —1, and we set 

Q = Te p £i T -1 £i , (73) 

where the matrix P parameterizes the massive modes. By the very definition of what it means to be a massive mode, 
P commutes with EJ in both cases. To do the integral over P to leading order in l/N we may put P equal to zero in 
the symmetry-breaking terms with parameters a 2 and £ (we remind the reader that both a 2 and C are of order l/N). 
Thus we need to integrate 

SDct Ar (Q)e- (Ar/2)STr(QSl)2 = (J N c NSTrP-(N/2)STrc 2p = a N Q - (N/2) -2N STr P 2 + . . . _ ^ 



We see that the fluctuations of the massive modes are of the order P ~ 1/y/N. Because of the smallness of these 
fluctuations we may replace the nonlinear Berezin measure for the P- variables by the flat Berezin measure DqP (i.e., 
the product of differentials for the commuting variables and derivatives for the anti-commuting variables). Thus we 
have DQ ~ DqP D^(Q a ) where Q a = TEjT _1 Ei and D^{Q a ) is the Berezin measure which is invariant under 
the transformation QEi 1— ► TQEiT -1 . The integral over the massive modes P then is a simple Gaussian integral 
/ DqP e~ 2N \ Doing it we immediately arrive at the result (f3"5|) of the previous section. Moreover, we are now in 
principle able to determine the precise normalization constant. We will insert the correct overall normalization when 
evaluating the partition function below. 
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V. EXACT CALCULATION USING SUPERBOSONIZATION 



We now use the result (|69p from superbosonization to derive an exact expression for finite N. To that end we start 
from the formula for the superdeterminant, 

SDetQ = SDetp^=^L, (75) 

where "qX -1 ^ means the scalar which is obtained by sandwiching the matrix X^ 1 between the row vector r\ and the 
column vector £ . In view of Eq. (|75p . the factor SDet 7V (Q) of the integrand of (|6"9")) is much simplified by making a 
shift y — > y + r/X' 1 ^ . Such a shift leaves the integral over y £ U(l) invariant: f(y) dy = f(y + Q dy . 
After this shift, our integrand depends on the anti-commuting variables only through the following factor: 

$ = c N,(A +{ i-a^v)x-^ + f(i-a-) {v x-^) 2 , A = A{X ) = ai X - ia 2 a 2 X - iz} . (76) 

Using the relation ^(^A -1 ^) 2 = Det _1 (A) ?7 + £ + ?y_£_ we now carry out the integral over the anti-commuting variables 
to obtain J d^+drj+d^-drj- $ = Det _1 (X) F(X,y) where 

F(X, y) =N(l-a 2 )+ N 2 (Det A(X) + (1 - a 2 ) y Tr A{X) + (1 - a 2 ) 2 y 2 ) . (77) 

We insert this into the integral representation (|6"9")1 of the partition function to get 

■-N-1 M N 

- I d 4 X Bct N - 2 (X) <j> 

J(i) 



Z^(z}\z,z*;a) = / d*X Det N ~ 2 (X) & dy y~ N+1 F(X,y) 

4tH(N- 1)! J x>0 /u(i) 

x g-eNTrX+iJVTr (Xcr 1 Rcz+Xa 2 Im z)-iNz}y- ^Tr {Xa 1 Xa 1 +a 2 Xa 2 Xa 2 ) + ^ (1-a 2 ) y 2 ^yg^ 

Next we do the y-integral. For this purpose we introduce the a-dependent functions 

h N , k (r) ■■= { dy{iyr N + 1 + k e-*v+*< 1 -W , (79) 

*K Al(l) 

which will be shown presently to be scaled Hermite polynomials. With this definition, what remains to be done is an 
integral over the positive Hermitian 2x2 matrices X : 

atN+2 r 

g _ IV / ^4y- e -eNTrX+iN'Tt{Xa 1 Rcz+Xa 2 Iniz)-^Tr(Xcr 1 Xcr 1 +a 2 Xa 2 Xa 2 ) 



2TT(N-iy.j x>Q 

x Bet N - 2 (X)([h Nt2 {Nz* f )+iTrA(X) h Nfl (Nz* f ) - (Bet A(X) + N' 1 ^ - a 2 )) h Nfi (Nz* f )\ . (80) 
To compute the A-integral one may use the parametrization 



x=f cU ^ + v2 + w2 u \ (u,v, w eR, P eR + ). (81) 

\ v + iw e \fp + v z +w z J 

The integration measure in these coordinates is expressed by 

d 4 X = 2dpdudvdw , (82) 

and some traces appearing in the exponent of the integrand are 

|Tr(X(7iXo-i) =p + 2v 2 , ±Tr(Xa 2 Xa 2 ) =p + 2w 2 . (83) 



A notable feature here is that the variable u occurs only in the factor e eArTrX = e 2tN Vp+ v2 + w2 cosh u . Thus the 
integral over u for fixed (3 :=2eN^/p + v 2 + w 2 ^ can be carried out and yields the hyperbolic Bessel function 

e- 0coshu du = K o (l3) . (84) 
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Let us now show how the functions ftjv 1 fc('?~) are expressed in terms of the Hermite polynomials H n (x) defined by 

H n {x) e~ x = (—1)™ d n /dx n e~ x . By shifting and setting the variable to zero after differentiation we can rewrite this 
definition as 



H n (x) = ( -^— e - y2+2x y 
\dy n 



= S^/ dyy- n - ic-^v+v*, (85) 

y=0 2m 7u(l) 



where we have used Cauchy's formula (d n f /dy n )(0) = (2iri) 1 n\ §un\ f(y) V n 1 dy . We ultimately want to take the 
limit N — > oo. To get a good view of the large- TV asymptotics we introduce the scaled Hermite polynomials 

h n (r) = (-1)" (2 c/n)""/ 2 nr'Hnir/V^n) = I dy y -n-l e -iTy+(n/2)(y 2 -l) (8g) 

2m Jv(i) 

By a saddle-point computation of the U(l) integral, these polynomials have the large-n behavior 

h n (r) ~ (nir)- 1/2 cos(r + mr/2) . (87) 
Comparing the integrals (|86[) and (|79|) we read off the relation 

n = N - 2 - k . (88) 



~h Ntk (r) = (1 - a 2 ) fc +"/ 2 (e7V/n)"/ 2 ft„ g2 



l2 ) 

A further simplification of Eq. (|80p is now achieved by the 3-term recursion formula 

Nh N , 2 (T) + rh NA ( T ) + (N - 2)(1 - a 2 )/tjv,o(T) = , (89) 

which results from partially integrating (N(l — a 2 ) y — d/dy) e^ 1- " = against c~ 1Ty y~ N+2 dy . Using the identity 
to eliminate the /ijv,2 term from (|80[) we arrive at 

Z_l = \7 / d PP N ~ 2 e- N{ - 1+a2 ^ / dv / dw e 2i7V(,Re Z +„I mZ )-2iV(, 2 +a^ U 2 ) ifo(2eArv / p+ ^2 + ^2) 



tt(JV-1)!7 

x ((2iu + 2a 2 w + 2 ;.)/i A r a (iVz}) + (z* f (2iv + 2a 2 w + z* f ) + p(l -a 4 ) - (1 -7V _1 )(1 - a 2 ))^v, (iVz})) . (90) 

This expression for the partition function is exact for all matrix dimensions N > 2 . (It is, however, false for N = 1 
because the superbosonization formula fails in that case; see the discussion in [2*ol, l2l|.) 



A. Calculation of the log e term 

We now extract from the integral representation (|90p the term which is singular in the limit e — > . For e — > we 
may replace the hyperbolic Bessel function by its leading logarithm, 

K (p) ~-logC9/2) = -loge-log (aVp + « 2 +™ 2 ) , (91) 

where we keep only the singular term loge for now. The integrals over the variables v and w then become Gaussian 
with mean values (v) = ^Rez and (w) — ^^jlmz and variances var(t>) = (4iV) _1 and v&v(w) — (4iVa 2 ) -1 . The 
remaining integral over p after scaling p — > A _1 (l + a 2 )~ 1 p yields the gamma function T(N — 1) = f °° p N ~ 2 e~ p dp 
and a similar term with N — 1 replaced by N. Altogether we obtain 

N 2 -f Rc 2 ( Z )-f Im 2 ( 2 /a) 

lim | log, \-'Z^ = {N _ 1){1 + a2)N - la (*} ~ z*){h Ntl (Nz}) + z) h N , (Nz})) . (92) 

In view of the asymptotic behavior ([87]) it is clear that this will tend to a good limit for N — ► oo when the product 
Nz*j is kept fixed. Inserting the definition of the polynomials h nt f. and using the recursion relation 

H n+1 (x) = 2xH n {x) - 2nH n - l {x) (93) 

with n = N — 2 , we find the simplified expression 

lim | loge T 1 ^! = a-'iz} - z*)e-^ Rc2 ^-^ lm ^ z /^ C N (a) (-if H N ^(bz}) , (94) 
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where the normalization constant Cjv(a) and scale factor b = b{a) are given by 

N N I N 

This closed-form expression for the log e contribution is still exact for all matrix dimensions N > 2 . 

Let us check that this result is consistent with the expression (|58[) obtained in the large- N limit. For that we 
observe that the Hermitc polynomials H n (x) for yfnx fixed and mco are asymptotic to 

(-l) n H n (x) ~ \/2(2 7Vc)" /2 cos(V2^a; + n7r/2) . (96) 

Recalling that in the microscopic limit we send A^ — > oo while keeping A^a 2 , A^z* and Nz*f fixed, we then see that the 
microscopic limit given in (|58[) is precisely reproduced. 



B. Contribution of order e 1 



We now set Zj — z* and compute the e° contribution to the partition function. This contribution is given by the 
log(p + v 2 + w ) term in the expansion (|91|) of the hyperbolic Bessel function. (From the preceding section we know 
that the constant terms in the expansion of Kq(P) yield zero for z*j = z* .) To facilitate the computation, we write 

- 2 log + v 2 + w 2 = lim (J ye- r( - p+v2+w ^+logS + ^, (97) 

where 7 is Euler's constant. The singular constant log S and 7 can be dropped as they make no contribution for 
= z* . The integrals over v, w, and p can then be carried out as before, and the resulting limit S — > exists. The 
order e° contribution is thus given by 

N -1 

Z_i(z*|z,z*;a) = (f(l-a 2 ))~ / dr 



(1 + a 2 + N-h-) N ~ 1 v / (2N + r)(2A^a 2 + r) 
N N 3 /N{ Rez ilmz \ (A r -l)(l-a 2 ) 



x lUte-^J^- W^ Fjv - 2M J ■ (98) 

This is the exact finite- N result for the bosonic partition function. Using the orthogonal polynomial approach it can 
also be obtained from the Cauchy transform of the fermionic partition function [331 ] . 

Let us take once again the microscopic limit (N — > 00, with Nz* . Nz and A^a 2 fixed). Doing so in Eq. (|9"8|) we get 

ijV-§JVa 2 roo ~(2Na 2 +r)- 1 Im 2 (JV z) /H m (M z ) \ 

Z-! = =— / dr OAr \ ' sm(Nz* + Nir/2) + cos(Nz* + JVtt/2) . (99) 

V2^ Jo e r V2Na 2 + r \2Na 2 + r y 1 J \ 1 , \ > 



To check this result we invoke the following identity: 

„-JV 2 lm 2 (z)/(2Na 2 +r) r Ar, 

dr 2 UA = / ^_ e _„r 8V +M, ta , (100) 

e'V2A^a 2 +r Jm 1 + q 2 

and a second identity of the same kind which is obtained by differentiating both sides of (jlOOp with respect to Im z . 
Using Euler's formula cos# + isin# = e lS we then immediately recover Eq. (|61| . 



VI. CALCULATION OF Z_i USING COMPLEX ORTHOGONAL POLYNOMIALS 

In this section we first derive the partition function with one bosonic quark from the Cauchy transform of orthogonal 
polynomials. In this approach no rcgularization procedure is required. For comparison with the er-model we also 
compute the partition function with one fermionic quark and two conjugate bosonic quarks, which diverges as loge . 

To apply the method of complex orthogonal polynomials to the FKS model we first express the Gaussian probability 
distribution for H and A given in Eq. @ in terms of the eigenvalues z^ of H + A and z^, of H — A . The joint 
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distribution for the eigenvalues was calculated in [12j. Including the expression for the exponent, which follows from 
the decomposition 

TtHH^ + ^TrAAt = 1 + 1 / a Tr (H + A)(H* + A+) + ^—^-Tv({H + A) 2 + (fft + A*) 2 ) , (101) 
the joint eigenvalue distribution function is given by fl2| 

P({z k ,z* k }) = C]l\z k - 2i | 2 e-T(i+i/« 2 )E fe W 2 -f (i-V» 2 )£*(4+4 2 ). (i 02 ) 
fc<; 

The partition functions (and correlation functions) of the FKS model can now be derived by means of the method of 
complex orthogonal polynomials with polynomials p n (z) defined through 

d 2 z u 

and weight function given by [4(| 

W (z,2*;a) = e-T( 1+1 /- 2 )l-l 2 -f(i-i/ Q 2 )(^+^ 2 ). (104 ) 
These polynomials, which have been known for some time (34l. l35l. [36j . are given by 



p»(z) ~ H n (bz), b =\j 2{l N _ a 2 ) ■ ( 105 ) 

where -ff„ are the Hermite polynomials. The p k are in monic normalization with respect to z . The leading a-depen- 
dence of rfe for large k is given by 

r k ~e^ ka ' 2 . (106) 

In this section, we will not keep track of numerical and a-dependent prefactors. 

General expressions for partition functions in terms of complex orthogonal polynomials have been given in [37l . l38l . 
[39j . Below we derive the explicit expressions for the microscopic limit of Z_i(z*\z, z*; a) and Z_\(z*f\z, z*\ a). 

A. The partition function Z-i(z*\z, z*; a) 

The partition function with one boson, Z-i(z*\z, z*; a), can be expressed as a Cauchy transform [37l. l38l, [39| 

Z^(z*\z,z*;a) = — !— [ d 2 z' w(z', z'*; aW-iO'*)— ^— ; . (107) 
rjv-i 7 z — z' 

In the microscopic limit where Nz and iVa 2 are kept fixed for N — > oo, the Hermite polynomials can be replaced by 
their asymptotic limit (]96[) and the weight function reduces to 

»Im 2 (i] 

u>(z,z*;a) ~ e . (108) 



We t hus find 



Z^^lz^V) ~ r-n^-ie-i^- / dx'dy' e"^- + 1 , J , f . (109) 

z — ar — iy' 



The integral over a;' can be performed by a contour integration, whereas the remaining integral over y' can be expressed 
in terms of the complementary error function. This leads to the expression 

Z^(z*\z,z*;a) ~ (-i) N eS Na ' f dy' e"^ (e iNz e 2Ny ' 6(-y' + Imz) + (-l) N e - iNz e- 2Ny ' 9(y' - Imz)) 



3^.-^(^/2) erfc (W-ATmzN e|jVa% _ iW(x+7r/2) erfc / 2iVa 2 + ATmz\ 



in agreement with (1631) . 
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B. The partition function Z-i(zf\z, z*\ a) 

In (l9j it was shown that the singular part of the chiral random matrix partition function with a pair of conjugate 
bosonic quarks and Nf fermionic flavors factorizcs. The same reasoning can be applied to the FKS model resulting in 

Z™{z}\z, z*;a)~ (z) - z*) z[ N ~V (z); a) z£l h Jz, z*;a) + 0( c °), (111) 

where the phase quenched bosonic partition function, Z^_ hos , is defined in j5|) and z[ N is the partition function 
with one fermionic flavor. Z^}_ hos is given by the weight function times log e 0, [l9[ 

Z pq - hos {z,Z ;a) = |loge| , (112) 

and the JVy = 1 theory can be expressed in terms of the orthogonal polynomials [l5[ 

z[ N - l \z* f] a)^p N ^{z* } ) . (113) 
Inserting this result and (|112[) into (|llip we obtain 

Z^\z}\z,z*\a) ~ (z}-z*)w(z,z*;a)e-i Na2 H N - 1 (bz* f ) \oge + O(e ) , (114) 

in agreement with the exact finite- ./V result in The factor exp(— 3Na 2 /2) results from a factor exp(— Na 2 ) from 
rjv-i and a factor exp(— Na 2 /2) from the ratio of pn-i (in monic normalization) and -ff/v_i . In the microscopic limit 
this results in 

Z- X (Z}\Z, Z*;a) ~ (z* f ~ Z*) q-^o? e -Nl m \z)/{2a 2 ) sin(A r z * + jy^/g) | loge| + O(e ), (115) 

in agreement with the result (|58[) obtained earlier in this paper. 



VII. CONCLUSIONS 



We have analyzed the bosonic partition function of a Hermitian random matrix model deformed by a nonhermitian 
random matrix model. We have shown that the microscopic limit of the partition function can be obtained by 
essentially using symmetry arguments only. There are, however, several subtleties that deserve attention. First, 
the partition function has to be regularized by multiplication with a conjugate bosonic and conjugate fermionic 
determinant. Second, because fermionic degrees of freedom are present, two incquivalent saddle-point manifolds have 
to be taken into account. Third, convergence of the partition function leads to the boson-boson block of the manifold 
of the Goldstone degrees of freedom being a noncompact subset of the set of positive definite matrices. 

The main advantage of the symmetry approach is that it gives a clear view at universality. Goldstone modes, 
which are separated from the rest of the excitation spectrum by a mass gap, decouple in the microscopic limit, and 
their mutual interactions are completely determined by the symmetries and the pattern of symmetry breaking of the 
microscopic partition function. This means that our results for the FKS model are valid for the microscopic limit of 
any model with the same symmetries and a mass gap. 

To obtain results for finite-size matrices one has to perform a detailed calculation. We have presented results 
using two different methods: the superbosonization method and the complex orthogonal polynomial method. The 
disadvantage of the orthogonal polynomial method is that universality is not manifest at all stages of the calculation. 
In the superbosonization method the universal partition function is obtained after integrating out the massive modes, 
which is a trivial step when the proper coordinates are used. The orthogonal polynomial approach, which is applicable 
to invariant random matrix models, has as its main advantage that it can be generalized in a straightforward way to 
any number of flavors. We have also performed the calculation using a hybrid method where the four-fermion term 
is decoupled by means of the Hubbard-Stratonovich transformation. Since this calculation did not provide additional 
insights, we have refrained from presenting it in this paper. 

For the present problem, the superbosonization approach does not have a clear advantage over the hybrid method, 
but in general we expect that it will be simpler to integrate out the massive modes if fermions and bosons are treated 
in a unified way. We also wish to stress that a major advantage of the superbosonization method is that it can 
deal with nongaussian probability distributions. Such distributions have important applications in, e.g., quantum 
gravity and growth phenomena. However, nongaussian perturbations do not affect the universal results obtained in 
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the microscopic limit. These are determined by symmetries and can be derived from symmetry arguments alone as 
we have shown in this paper. 
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